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Abstract 

vol 

£\\ , The resonances (r]c,XcO,c2) contribution to two-photon exchange (TPE) effects in e + e <-» pp is 

calculated in a simple hadronic model. The calculation shows the TPE contributions by resonance 
Xc2, which are dependent on the unknown phases 4>e,m of proton's time-like form factors Ge,m, are 
much larger than the TPE contributions by non-resonance and are comparable with measurement 
precision of coming PANDA detector at y/s ~ M x 2 for most <j>E m- 

^ . PACS numbers: 13.40.Gp, 13.60.-r, 25.30.-c. 
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1 Introduction 

£NJ 1 The structure of nucleon is one of the most important topics in the hadronic physics. The electromagnetic 

interaction provides a clean method to measure such a structure. One of the most important measurement 
is to determine the electromagnetic form factors Ge(Q 2 ) and Gm(Q 2 ) of nucleon. The measurements of 
R = fj,pG e(Q 2 ) I 'G m(Q 2 ) at space-like region by polarized method and Roesenbluth method [T] indicate 
it is not trivial to extract the form factors directly from the angle dependence of cross section. Since then, 
the two-photon exchange (TPE) effects attract many interest. In the literature, the TPE contributions 
in ep — > ep and e + e~ <-> pp have been estimated by some model dependent methods pHl] and model 
independent analysis [8HT5] (see recent review articles [2]). Experimentally, to detect the TPE effects 
directly, the measurements of R e + e - defined as the ratio of e + p to e~p differential cross sections are 
proposed by VEPP-3 [15], JLab [16] and OLYMPUS [17] . and the measurement of forward-backward 
asymmetry in pp —¥ e + e~ is proposed by PANDA [18J. In this work, we estimate the TPE contributions 
in the unpolarized processes e + e~ O pp when y/s lies in the resonances T) c , Xc0,c2 region. We arrange our 
discussion as follows: in Section 2, we review the contributions of vector resonances (ip, in e + e~ — > pp 
by one-photon exchange (OPE); in Section 3, we discuss the TPE contributions by ?y c ,XcO,c2 and in 
Section 4 we present the numerical results and give a discussion. 
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2 Resonances contribution in e + e — > pp by OPE 

Considering the process e + (/c2) + e _ (fci) — » p(p2) + p(pi) at tree level, by the vector-meson-dominance 
model, the 1 resonances contributions can be described by diagrams in the left panel of FigQ] where 
we limit our discussion in the charmonium region. The couplings of 1 charmonium ip,^' with pp are 
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Figure 1: Left panel: diagrams for 1 resonances (ip,ip') contributions in e + e _ — > pp by one photon 
annihilation; right panel: diagrams for h ,0 ++ ,2 ++ resonances (?7c,XcO,e2) contributions in e + e~ — > pp 
by two photon annihilating. 



described as in Ref. [19] with a little different notation 



rtoM'rf) = - l[gM + ^gW), (1) 

where Mjy = 0.938 GeV is the mass of proton, q = (y/s, 0,0,0) in the cm. frame [20]. By these Born 
diagrams, the unpolarized differential cross section in the cm. frame can be written as 

d*°_ a 2 em M*(l - 4M 2 f / S ) 1 / 2 [AMl\G%\ 2 (\ - cos 2 g) + gjG&Pft + cos 2 0)] 

dn 4,2 /|((«-M|)2 + r|M|) 

«L^(1 - 4M 2 f / s ) 1 / 2 [4M& |G%'| 2 (1 - cos 2 0) + .s|Gt;| 2 (l + cos 2 0)] 
4* 2 / 2 ,(( s -Af 2 ,) 2 +r 2 ,M 2 ,) 

a 2 m ^M 2 ,(l - 4^/^ 4M*G|g^ - cos 2 g) + sC^G^l + cos 2 g) 

4s 2 61 Uf^is-M^ + iT^M^is-M^, -tiyAfy) h U 



with M.0 = 3.096 GeV, Ml = 3.686 GeV the masses of V>, V^i $ the scattering angle in the cm. frame, and 
it is a good apprc 

by QI] 



G^'^ = 9i>,i>' + K i/j,i/j's/4Af^r, Gjy- = g^.if,' + K^^/. Since we limit the discussion in the charmonium 
region, it is a good approximation to treat the couplings G^^ as constants which are constrained 



fl _ AM 2 IM 2 W 2 
T < ,, - - * m Nl m i>,psi') (9M 2 , M 2 |/7<M>'|2\ 

_teL« M _ l-4My s |G^'/^'| 2 ... 

with Otyrf ~ 0.67 at corresponding y/s = M^^t. In the following calculation, for simplicity, we naively 

assume that the phases of G^ E and G^ , GfjL and are the same, respectively. By these assumptions and 

the constraint conditions, only two physical parameters 4>m,e which are defined as G%* u = \G%* U \ 
are left as unknown. 



3 Resonances contribution in e + e — > pp by TPE 

When considering the TPE contributions from other resonances in charmonium region, the diagrams 
shown in the right panel of Fig Q] should be considered, where only h , ++ , 2 ++ resonances are included 
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since their couplings to pp and 2j are relatively large. The couplings of + , ++ , 2 ++ charmonium with 
two-photon can be described as in |21) . 

= -i9^ 11 De^^{k 2 + q)^(k 1 -q) V21 
r^ oo = -* 5xc077 # 2 {[(fci - ■ (fe + q)9^ - ~ q r (fe + gH^ + (^1 " 9) • (fe + g)] 

+ (fc 2 + q)^ (k 2 + (fcl - q? + (*i - (fci - g)" 1 (fc 2 + <z) 2 

- (fci - <z) 2 (fc 2 + g) V m - (fc 2 + q) Vl (fci - g)" 1 (fci - g) • (fc 2 + g)}, 
r 77™ = -^x.77^ 2 {/ lfI2 »^ ! (fci - g) • (*a + g) - a" 1 " 1 (fci - g)" 2 (*a + g)" 2 

- (fc 2 + g)" 1 (fci - g)" V 1 " 3 - (fci - g)" 1 (*a + g^V 1 " 2 }, (4) 

where k\ — q,k 2 + q are the momentums of photons, ^ii,v\ are the corresponding Lorentz indexes and 
D = [(q + fc 2 /2 — fci/2) 2 — m 2 + ie] _1 with m c the mass of c quark, and we take m c — 1.5 GeV in the 
following calculation. 

For the couplings of ^,0 ++ ,2 ++ charmonium with pp, for simplicity we take them the same as 
those with e + e _ , which can be described as in |21j . 

?)cPP 9r) c ppl5 > 

T x"c2pp = 9xc2 P p(Pi -ViYl*- (5) 
These couplings are constrained by the physical decay width T(ij c ,XcO,c2 — > 77, PP) with 

,2 M a 



r 



?77Xc2 Xc2 



X=2-s-77 



lS^m^ + M 2 ^' 



9l* P pJ M L ~ 4M %( 3M L ~ * M L M n ~ 32M 4 



T Xc^pp ~ 120ttM2 c2 ' (6) 

where M X2 = 3.556 GeV is the mass of x c2 , r Xc2 ^ 77 = 1.97 MeV x(2.56 x 1CT 4 ) = 5.04 x 1(T 4 MeV 

"< _ — 1 07 A/TWU v f7 O v m-S^l _ 1 /I' 

Xc2-*-PP 



and r Xc2 ^ P p = 1.97 MeV x(7.2 x 1CT 5 ) = 1.42 x 1CT 4 MeV [22]. This results in 



IS77Xc 2 1 = 1-57 x lCr 2 , \g Xc2PP \ = 7.39 x 1(T 4 . (7) 

Here, only g Xa2P p and <? 77Xc2 are given since the calculation shows the contributions from r\ c and XcO are 
identically zero because of the Dirac structure of corresponding unpolarized differential cross sections. 
For simplicity, we assume the phases of these two couplings are zero. 

The propagator of \ci is described as the standard Breit-Wigner form [23] 

A 1 pviP^ DM2 pp pv2 pp pH2 p" 

Qfl2f2;pU _ i — ln^Pt I „VlU I „P2U I „^2P \ 

x "2 p2 _ np +iM T 2 M 2 M 2 M 2 M' 2 

1 O 1 PM2 p"2 1 pp DW 

+ 2^ P 9 V2U + 9^9 U2P g^gn + -(-g^ - -^-){-gP" - ^)}- (8) 



With Eqs. (|4I5I8[) . the corresponding amplitudes for the diagrams in the right panel of FigfT] can 
be written down. Their interferences with Born diagrams can be calculated directly and we use the 
FeynCalc [23] to do the analysis calculation and LoopTools [25] for the numerical calculation. 
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Figure 2: The cos6* dependence of 8 at y/s — M Xc2 by \c2 contributions to TPE effects. The black solid, 
red dashed and green dotted curves are corresponding to <j)E — 0, 7r/3 and 2/3tt, the blue dashed-dotted 
curves named as EPJA are explained in the text. 

4 Results 



Since only Xc2 gives the contributions, we limit the discussion in y/s ~ M Xc2 . To show the TPE contri- 
butions, we define 

da^ /dn 

A[go _ gi] _ 2tt (da(°) + da^)/dQsin9d9 - 2irf*J°(daW + da^)/dQsin9d9 
2TrJ^(daW+daW)/dnsin9de + 27rj^g"(dcTW + daW) / dtlsinOdO 
Jg 1 daM /dnsinOdO 



(9) 



with da^/dQ, the unpolarized differential cross section from the interference of TPE diagrams and 
Born diagrams. By the crossing symmetry and Lorentz invariance, we have the relation for the am- 
plitudes M(e+(fca)e-(fci) -> p(pa)p(pi)) = M(p(-pi)p(-j»2) -> e + (-*i)e-(-*a)) = M{p{ Pl )p(p 2 ) -> 
e + (fci)e _ (/c2)), which means the ratios of cross sections in e + e — > are the same as those 

in pp — > e + e~ . 

The 9 and </>m dependence of 8 at y^s = M Xc2 are presented in Fig[2]and FigEJ The results show the 
TPE contributions to 8 by \c2 are almost independent on <f>E, especially at small 9, while are strongly 
dependent on 9 and 4>m- The corrections 8 are odd functions of cos 9, which is a general property of 27 
annihilation effects. The most interesting property of Figj3] is that the absolute corrections \8\ are larger 
than 10% at small 9 (9 < tv/10) for almost all <f>M- 

To compare with the measurement precision of coming PANDA detector, we present the corresponding 
8 used by the PANDA simulations [35] as the curves labeled EPJA in FigJ^l By the definition of our 8 
and TPE contribution Eq.(14) of Ref. [25], w e have 



2^t{t + 1){G e /t - G M )F 3 cos(9)sin 2 {9) 



GUl + co S 2(9)) + (G 2 E /T)sin 2 (9) ' (10) 



f2 



with r = s/AM 2 N and F 3 defined by Eq.(14) of Ref. [26]. The curve of 8 named as EPJA in Fig [2] is 
corresponding to the case F 3 /Gm — 5% assuming Ge = Gm- The results in Fig J5] show that the realistic 
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Figure 3: The 4>m dependence of 5 at y/s = M Xc2 by \c2 contributions to TPE effects. The black solid, 
red dashed, green dotted and blue dashed-dotted curves are corresponding to 4>e = 0, 7r/3, 2/3tt and it. 



8 in the small 9 {\cos9\ > 0.8) are much larger than 5(EPJA) for most 4>m- In the region where PANDA 
detector works, the realistic S are larger than (5(EPJA) in the region 0.5 < \cos6\ < 0.8 and comparable 
with J(EPJA) in the region \cos9\ < 0.5 for some 4>m, for example 4>m = 1/3tt, 2/3tt. 

The full results show that the TPE contributions to 5 from \c2 are much larger than the usual TPE 
effects from non-resonance contributions. The latter are usually less than 1% at small 9 [51[7] and is a 
challenge to be observed by PANDA [7] , while the former are comparable with the measurement precision 
of PANDA at yfs ~ M X c2 • This suggests that the coming PANDA experiment may detect the direct TPE 
effects at y/s ~ M Xc2 ±0.0f GeV. And outside of this region, the TPE contributions from \2 are as small 
as 1% at the small 9 (9 < n/10) and the corrections to the measurements of Gm,e will be very small [26] . 

Fig |4]displays the curves of forward-backward asymmetry A^^ 6 ^ vs. 4>m at ^fs — M Xc2 , M Xc2 ±0.001 
GeV with 4>e = as an example. The left panel of Figf?] shows that the forward-backward asymmetry 
with 9q — 7r/5,#i = 7r/2 is larger than 2% for almost all <f>M when combining the three y/s points. 
The right panel of Fig f4] shows that the forward-backward asymmetry with smaller 9\ = 7r/3 is much 
enhanced. Fig (5] displays the 4>e,m dependence of ^I^/ 5- ^/ 3 ] and J 4l 7r / 5_7r / 2 ] in contour form, which shows 
the properties of forward-backward asymmetry in the parameter space more clearly. By the simulation 
of PANDA detector [25] , the reconstructed number of counts at s = 12.9 GeV 2 (~ ^- Xa2 ) is about 10 3 for 
\cos9\ < 0.8 (corresponding to 9 > l/57r). The left panel of FigJS] shows that the number of asymmetry 
events for 9 C [ir/5 — ir/2] at y/s ~ M Xc2 is larger than 20 for almost all 4>e,m and larger than 40 most 

(f>E,M- 

To show the property of asymmetry events number, we define the ratio R as 

\l x do-W /dnsin9d9 

R=^- 2 . (11) 

fe f daW /dttsin6d6 

The angle 9\ dependence of ratio R at y/s = M Xa2 with = tt/5 is presented in FigJS] which shows two 
interesting properties: (1) R is almost independent on 4>e,m except for at very small 4>m- The special 
behavior of R at very small 4>m is corresponding to the property of 5 at </>m = in Figj2j which is not 
monotonic function of cos9 in the region 9 C [tt/5,tt/2]; (2) R reaches about 80% for 9\ = tt/3 and 
90% for 9\ = 2tt/5. These two properties indicate that whatever 4>e,m are, for experiments, most of the 
forward-backward asymmetry events lie in 9 C [n/5, 7r/3]/[27r/3, 4tt/5] when the detector works in the 
region \cos9\ < 0.8. 

In conclusion, our calculation shows the resonance \c2 contributions to TPE effects may give large 
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Figure 4: The <j>M dependence of forward-backward asymmetry at y/s ~ M Xc2 with 4>e — 0. The red 
dashed, black solid and blue dotted curves are corresponding to y/s — M Xc2 + 0.001 GeV, M Xc2 and 
M Xc2 - 0.001 GeV. 




Figure 5: The 4>m.e dependence of absolute forward-backward asymmetry at y/s = M Xc2 in contour form. 

corrections in e + e~ O pp when y/s ~ M Xa2 . The results suggest that the coming experiment PANDA 
may detect the direct TPE effects in small scattering angle region such as 9 C [ir/5,Tr/3] at y/s ~ M Xe3 . 
For lower energies around 2.3 GeV, the contributions from resonances of light flavor quarks should also 
play similar important role in e + e~ ■(-> pp. 
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Figure 6: The 9\ dependence of R at y/s = M Xc2 with 9 = ir/5. The black solid, red dashed, green 
dotted and blue dashed-dotted curves are corresponding to 4*e — 0, tt/3, 2/3tt and n. 

References 

[1] R.C. Walker et al., Phys. Rev. D 49 (1994) 5671; L. Andivahis et al, Phys. Rev. D 50 (1994) 5491; 
M.K. Jones et al, Phys. Rev. Lett. 84 (2000) 1398; O. Gayou et al., Phys. Rev. Lett. 88 (2002) 
092301. 

[2] P.G. Blunden, W. Melnitchouk, J.A. Tjon, Phys. Rev. Lett. 91 (2003) 142304; S. Kondratyuk, P.G. 

Blunden, W. Melnitchouk, J.A. Tjon, Phys. Rev. Lett. 95 (2005) 172503. 
[3] Y.C. Chen, A.V. Afanasev, S.J. Brodsky, C.E. Carlson, M. Vanderhaeghen, Phys. Rev. Lett93 (2004) 

122301. 

[4] Dmitry Borisyuk, Alexander Kobushkin, Phys. Rev. C 78 (2008) 025208. 

[5] Dmitry Borisyuk, Alexander Kobushkin, Phys. Rev. D 79 (2009) 034001; Nikolai Kivel, Marc Van- 
derhaeghen, Phys. Rev. Lett. 103 (2009) 092004. 

[6] D.Y. Chen, H.Q. Zhou, Y.B. Dong, Phys. Rev. C 78 (2008) 045208; Hai Qing Zhou, Dian Yong 
Chen, Yu Bing Dong, Phys. Lett. B 675 (2009) 305. 

[7] Julia Guttmann, Nikolai Kivel, Marc Vanderhaeghen, Phys. Rev. D 83 (2011) 094021. 

[8] Pierre A.M. Guichon, M. Vanderhaeghen, Phys. Rev. Lett.91 (2003) 142303; 

[9] Michail P. Rekalo, Egle Tomasi-Gustafsson, Eur. Phys. J. A22 (2004) 331; E. Tomasi-Gustafsson, 
G.I. Gakh, Phys. Rev. C 72 (2005) 015209. 

[10] J. Arrington, Phys. Rev. C 69 (2004) 032201; J. Arrington, Phys. Rev. C 71(2005)015202; J. Ar- 
lington, W. Melnitchouk, J. A. Tjon, Phys. Rev. C 76 (2007) 035205. 

[11] Yu-Chun Chen, Chung- Wen Kao, Shin-Nan Yang, Phys. Lett. B 652 (2007) 269. 

[12] M. A. Belushkin, H.-W. Hammer, U.-G. Meissner, Phys. Lett. B 658 (2008) 138. 

[13] G.I. Gakh, E. Tomasi-Gustafsson, Nucl. Phys. A 771 (2006) 169; E. Tomasi-Gustafsson, E.A. Kuraev, 
S. Bakmaev, S. Pacetti, Phys. Lett. B 659 (2008) 197. 

[14] Carl E. Carlson, Marc Vanderhaeghen, Ann. Rev. Nucl. Part. Sci. 57 (2007) 171; J. Arrington, P. 
G. Blunden, W. Melnitchouk, Prog. Part. Nucl. Phys. 66 (2011) 782. 

[15] J. Arrington, D. M. Nikolenko et al., |nucl-ex/040 8020; D.M. Nikolenko et al., Phys. Atom. Nucl. 73 
(2010) 1322; D.M. Nikolenko et al., PoS ICHEP 2010 (2010) 164. 

[16] W. Brooks et al., Jefferson lab experiment E07-005. 

[17] M. Kohl, AIP Conf. Proc. 1160 (2009) 19; M. Kohl, AIP Conf.Proc.l374(2011)527. 

[18] PANDA Collaboration, arXiv: 0903.3905; M. Sudol et al., Eur. Phys. J. A 44 (2010) 373. 

[19] T. Barnes, X. Li, W. Roberts, Phys. Rev. D 77 (2008) 056001. 

[20] Where q is corrected as (i/s, 0, 0, 0) in cm frame comparing with [19] . 



7 



[21] Johann H. Kuhn, Jean Kaplan, El Ghali Oudrhiri Safiani, Nucl. Phys. B 157 (1979) 125; D. Ebert, 

R.N. Faustov, V. O. Galkin, Mod. Phys. Lett. A 18 (2003) 601. 
[22] K Nakamura et al (Particle Data Group)., J. Phys. G 37 (2010) 075021. 

[23] Shi-Zhong Huang, Peng-Fei Zhang, Tu-Nan Ruan, Yu-Can Zhu, Zhi-Peng Zheng, Eur. Phys. J. C 
42 (2005) 375. 

[24] R. Mertig, M. Bohm, A. Dcnncr, Comput. Phys. Commun. 64 (1991) 345. 
[25] T. Hahn, M. Perez- Victoria, Comput. Phys. Commun. 118 (1999) 153. 
[26] M. Sudol ct al., Eur. Phys. J. A 44 (2010) 373. 



8 



